Measuring eccentricity in binary black-hole initial data 
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Initial data for evolving black-hole binaries can be constructed via many techniques, and can 
represent a wide range of physical scenarios. However, because of the way that different schemes 
parameterize the physical aspects of a configuration, it is not alway clear what a given set of 
initial data actually represents. This is especially important for quasiequilibrium data constructed 
using the conformal thin-sandwich approach. Most initial-data studies have focused on identifying 
data sets that represent binaries in quasi-circular orbits. In this paper, we consider initial-data sets 
representing equal-mass black holes binaries in eccentric orbits. We will show that effective-potential 
techniques can be used to calibrate initial data for black-hole binaries in eccentric orbits. We will 
also examine several different approaches, including post-Newtonian diagnostics, for measuring the 
eccentricity of an orbit. Finally, we propose the use of the "Komar-mass difference" as a useful, 
invariant means of parameterizing the eccentricity of relativistic orbits. 

PACS numbers: 04.20.-q, 04.25. Dm, 04.70.Bw, 97.80.-d 



I. INTRODUCTION 



The possible detection of gravitational waves by detec- 
tors such as LIGO and LISA is driving rapid progress in 
the binary black hole (BBH) problem. The final stages of 
the inspiral and coalescence are believed to be primary 
sources of gravitational waves at frequencies accessible 
by such detectors. Theoretical models that accurately 
predict these final stages of inspiral are needed to help 
analyze and improve the rate of detection with future 
data. There are two techniques commonly used to study 
such systems, the post Newtonian approximation and nu- 
merical relativity. Since post Newtonian methods are 
expected to break down at smaller separations, accurate 
models of final stage inspirals need to be handled with 
numerical techniques. 

Numerical relativity breaks the problem into two 
parts, the specification of initial data and the numeri- 
cal time evolutions of this data. With recent advances in 
evolutions [H [3, S 0], it is as important as ever to fully 
understand the initial data one is starting with. Post 
Newtonian methods have shown that BBHs starting with 
large separation will evolve toward an adiabatic inspiral 
that follows a series of quasi-circular orbits Q. This qua- 
sicircular inspiral is thought to last till late stages where 
the black holes plunge toward coalescence. Most work 
to date in both the construction of BBH initial data and 
evolutions have focused on quasi-circular configurations. 
But the study of binaries in truly eccentric orbits near 
coalescence may be important for gravity-wave detectors 
(especially LISA) and, in any case, is of considerable the- 
oretical interest. A primary goal of this paper is to take a 
preliminary look at how current quasiequilibrium meth- 
ods for constructing BBH initial data can be extended to 
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construct general eccentric configurations. 

To predict quasicircular orbits, two techniques have 
commonly been used: the effective-potential (EP) 
method and the Komar-mass ansatz [3]. The Komar- 
mass ansatz compares two definitions of energy: the 
ADM energy[S] and the Komar mass[9]. The ADM en- 
ergy gives a proper definition of energy at spatial infinity 
in all cases, but the Komar mass is only accurate when 
a configuration is at least momentarily stationary. So, 
the Komar-mass ansatz posits that when these two ener- 
gies are equal, the system should be in quasiequilibrium. 
The EP technique is motivated by the Newtonian effec- 
tive one-body problem, and determines a quasi-circular 
configuration by finding models that have a minimum 
in the binding energy along a sequence of constant an- 
gular momentum, constant masses, and constant spins. 
The two methods have been contrasted and largely agree, 
with significant differences occurr ing only at separations 
close to final plunge of the BBHsflO. [ll|. 

It has been conjectured [Tot that the EP method, when 
used with methods for constructing quasiequilibrium ini- 
tial data, should generalize to represent BBHs in eccen- 
tric orbits. For initial-data methods designed to produce 
quasiequilibrium data, the generalization is not obvious 
because the individual black holes will not follow the in- 
tegral curves of any approximate Killing field. We will 
explore this conjecture extensively in this paper. In or- 
der for this conjecture to be true, the initial data used to 
construct the effective potentials must represent BBHs at 
turning points in their orbits. Unfortunately, we cannot 
test this property of the data directly. We will justify 
the assumption that the data represent BBHs at turning 
points and test it by comparison to various results from 
post-Newtonian methods. 

By constructing effective-potentials that represent 
BBHs in eccentric orbits, it becomes possible to make 
direct estimates of the eccentricity of an orbital config- 
uration. Of course, it is not possible to justify a unique 
definition for eccentricity for relativistic configurations 
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where elliptic orbits do not exist. Instead, we compare 
several different definitions [l^l and show that they give 
reasonable and consistent results. Because there is no 
unique definition of eccentricity, it is perhaps not the best 
quantity for parameterizing eccentric orbits. We find that 
the "Komar-mass difference" , the quantity that is set to 
zero in the Komar-mass ansatz for defining circular or- 
bits, may serve as a useful parameterization of eccentric 
orbits. 

We begin in Sec. HIl with an overview of the quasiequi- 
librium method we use to construct BBH initial data, 
emphasizing the aspects that will be most relevant to our 
subsequent discussion. In Sec. lIIH we will discuss various 
aspects of effective potentials in the context of their use 
with BBH initial data. In particular, we will justify as 
fully as we can the extension of these effective potentials 
to BBHs in eccentric orbits. In Sec. IIVI we will explore 
and compare several definitions of eccentricity, and will 
motivate the use of the "Komar-mass difference" as an 
invariant means of parameterizing BBHs in eccentric or- 
bits. All of the preceding discussions have dealt with 
non-spinning black holes. In Sec. |Vl we will briefly re- 
visit the case of corotating black holes. We end the paper 
in Sec. IVII with some conclusions. 



II. INITIAL DATA 

The binary black-hole initial-data sets that are used 
below were described in detail in Refs. [13, [3] and refer- 
ences within. Here, we give an overview of the methods 
used to construct the initial data, elaborating only on 
the details most relevant to finding black-hole binaries in 
quasi-circular orbits. 

Our initial-data sets are constructed within the ex- 
tended conformal thin-sandwich approach [isj. This 
approach is based on the standard 3 -I- I decomposition 
where the space-time interval is written as 

ds2 ^ -a^dt^ -f 7i^ (da;* + I3'dt){dx^ -f /3Mt). (1) 

Here, 7y is the spatial metric, and a and /3* are the 
lapse function and shift vector. Minimal initial data for 
a Cauchy evolution requires that we fully specify 7^ and 
the extrinsic curvature Kij (essentially a first time deriva- 
tive of the spatial metric) defined as 

Kij — ——Cnjij, (2) 

where n'' is the time-like unit normal to the t — const, 
initial-data surface. 

The conformal thin-sandwich (CTS) approach requires 
a conformal decomposition of jij , and that we specify the 
conformally related metric 7^ . In this work we will al- 
ways take 7ij to be flat. We must also specify the time 
derivative of the conformal metric dt"fij and the trace of 
the extrinsic curvature. We fix both quantities to be zero. 
The CTS approach then requires that we determine the 



conformal factor ^p relating 7^^ and the conformal metric, 
and the shift vector /3\ These are obtained by solving 
elliptic versions of the Hamiltonian and momentum con- 
straint equations (see Ref. [To] for details). Finally, the 
extended CTS approach also requires that we determine 
the lapse function a by fixing the time derivative of the 
trace of the extrinsic curvature and then solving the evo- 
lution equation for the trace of the extrinsic curvature 
as an elliptic equation. We fix the trace of the extrinsic 
curvature to be constant in time. 

In constructing black-hole initial data, we excise the 
black-hole interior from the computational domain and 
must impose boundary conditions at the excision sur- 
faces when solving the elliptic equations for ip^ and 
a. We demand that each black hole be in quasiequilib- 
rium by imposing the boundary conditions worked out in 
Refs. [l6| and [l3|. The assumptions of quasiequilibrium 
are essentially the same as those required of an "isolated 
horizon" (see Refs. [13, EM [3). We must also impose 
boundary conditions at the outer boundary of the com- 
putational domain (either at infinity or some large ra- 
dial distance from the black holes). For this, we assume 
that our configuration is asymptotically flat. However, 
asymptotic flatness does not fully fix the boundary condi- 
tions on the constrained data. The asymptotic condition 
on the shift is 

(3%^^ = (^0 X ry, (3) 

where JIq is an angular velocity vector. 

The time coordinate threading through our initial-data 
slice is defined by 

i^ = an''+/3^. (4) 

Imposing Eq. ([3]) , we see that fig determines the rotation 
of the "helical" time coordinate. For a binary system in 
equilibrium, the time coordinate generates a symmetry 
and the bodies move in circular orbits along integral lines 
of the time coordinate. For relativistic systems, the bi- 
nary can only be in quasiequilibrium and the time coordi- 
nate generates an approximate symmetry. Nevertheless, 
r2o represents the orbital angular velocity of the binary as 
measured by observers at infinity. From a computational 
perspective, ^Iq must be chosen. Ultimately, it is the ef- 
fect of different choices for JIq that we will be exploring 
in this paper. 

In Ref. [l3| , two independent methods for choosing the 
magnitude of flo were compared. Both methods attempt 
to produce a binary system that is in quasiequilibrium 
with the black holes in quasi-circular orbits. One method 
is based on the Komar-mass ansatz, first proposed by 
Gourgoulhon et. al.|2j, which posits that if is chosen 
so that the ADM energy J^adm and the Komar mass 
of a system are equal, then the system will be nearly 
stationary (i.e. in quasiequilibrium) and the binary will 
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be in a quasi-circular orbit. These masses are defined via 

E.u^ = / V, (gf ~ sfg) <fs\ (5) 

Mk = / (V.a - P'K,,) d?S\ (6) 

where Qij = jij — fij, fij is the flat metric, and Vj is 
the covariant derivative compatible with . The second 
method assumes that quasi-circular orbits are found at 
the minima of a reduced two-body effective potential. 
Effective potentials will be described in more detail in 

secnni 
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FIG. 1: Effective-potential plot for non-spinning equal-mass 
black holes constructed from numerical initial data 10]. The 
thin solid(multi-color) lines, are individual EP curves. Some 
of these EP curves are labeled it value of J/pm. Passing 
through the local minima of these EP curves and drawn as 
a bold(red) line is the EP sequence of quasi-circular orbits. 
The Komar sequence of quasi-circular orbits is displayed as a 
dashed(blue) line. 

In Ref. [l^l , we showed that the circular-orbit configu- 
rations produced by both methods agree remarkably well 
for both non-spinning and corotating black-hole binaries. 
Here, we simply show a figure that directly compares 
sets of circular-orbit models as we vary the binary sep- 
aration. Figure [1] displays the effective potential (EP) 
curves for non-spinning equal-mass black-hole binaries as 
computed in Ref. ^^T^. The vertical axes display the di- 
mensionless binding energy of the binary E^/ fi, where E^, 
is defined in Eq. ([S]), /i = mim2/rn is the reduced mass 
of the system and m = ml + m2 is the total mass of 
the binary. The horizontal axes display the dimension- 
less proper separation of the binary £/m. It is impor- 
tant to note that each EP curve consists of a sequence 



of models where the value of flo changes monotonically. 
Passing through the minima of the EP curves is the "EP 
sequence" defined as the sequence of quasi-circular or- 
bit models where quasi-circular orbits are defined via 
the effective-potential method. Also shown in Fig. [T] is 
the "Komar sequence" defined as the sequence of quasi- 
circular orbit models defined via the Komar-mass ansatz. 
It is clear that the two sequences agree quite well except 
for the smallest separations. See Ref. for a more 
detailed comparison. Most important for our consider- 
ations is that both the Komar-mass and EP methods 
choose particular models on each EP curve, with very 
similar values of flo, as quasi-circular orbit models. 

These sequences of quasi-circular orbit models repre- 
sent a good approximation to the adiabatic inspiral of a 
black-hole binary so long as the binary separation is not 
too small. For large enough separation, the time scale 
for radiation reaction to induce a significant change in 
the orbital radius is much larger than the orbital period. 
At small enough separations, radiation reaction becomes 
significant and the sequences of quasi-circular obits will 
eventually become a poor approximation to an adiabatic 
inspiral. It is difficult to estimate exactly where this tran- 
sition occurs, but it will certainly occur at separations 
larger than the inner- most stable circular orbit (ISCO). 
For a Komar sequence, the ISCO is defined to occur at 
the minimum in the binding energy. For an EP sequence, 
the ISCO is defined to occur at the point where the EP 
curves no longer have a local minimum. Both ISCOs 
occur near £/m ~ 5 (see Fig. [ij. 

When the CTS equations and excision boundary condi- 
tions as described above are used to construct initial data 
for binaries in quasi-circular orbit, the resulting data are 
consistent with a system in quasiequilibrium. In partic- 
ular, we have both bodies in the binary following the 
integral curves of an approximate helical Killing vector 
(the time coordinate). Furthermore, half of the initial 
data that must be specified is also consistent with this 
notion of quasiequilibrium. Recall that the conformal 
metric 7ij , the trace of the extrinsic curvature, and their 
time derivatives must be specified, and that we take all 
of the time derivatives to vanish. As long as the binary 
is in a circular orbit, the notion of quasiequilibrium is 
satisfied. 

However, we are left to wonder, what happens if we 
choose D,Q so that the binary is not in a quasi-circular or- 
bit? In Ref. [10|, we made the assertion that the resulting 
initial data would represent a binary at either pericenter 
or apocenter of a general bound or unbound orbit. To 
arrive at this conclusion, we are forced to give up the 
notion that the helical time vector represents an approx- 
imate Killing vector of the space-time. Our goal below is 
to see if this interpretation is reasonable. 
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III. EFFECTIVE POTENTIALS 



In Newtonian physics, 1-D effective potentials nicely 
capture the important features of certain dynamical sys- 
tems. In the case of the reduced gravitational 2-body 
problem, for a given orbital angular momentum, the ef- 
fective potential can be used to locate the turning points 
for an elliptic orbit of given energy, or the radius and en- 
ergy of a circular orbit. No such 1-D effective potential 
can be rigorously and uniquely derived for the fully rela- 
tivistic gravitational 2-body problem. However, a useful 
effective potential has been defined 0, [111, [lO, [IH in di- 
rect analogy with the Newtonian gravitational effective 
potential. 

In essence, an EP curve is the total energy of the sys- 
tem for a sequence of configurations where the radial sep- 
aration varies, but all other physical parameters are held 
fixed. In order to correspond to an effective potential, 
the velocity of the generalized coordinate that is allowed 
to vary (the radial separation in this case) must vanish so 
that there is no associated kinetic energy. Binary systems 
in bound orbits with vanishing radial velocity are either 
at apocenter or pericenter, and collectively we refer to 
these as turning points in the orbit. 

For black-hole binaries, we use the following func- 
tional definition for an EP curve. We consider a se- 
quence of initial-data configurations that hold constant 
the apparent-horizon masses of the individual black holes 
(toi and 7712)1 the magnitude and direction of the spins of 
the individual black holes, and the total angular momen- 
tum of the system"'^. We also require that the binary be 
at a turning point in the orbit. For the value of the effec- 
tive potential, we could simply use the total ADM energy 
of the system. However, it is more intuitive to use the 
binding energy since this allows us to estimate if a given 
configuration is in a bound orbit. The binding energy 
requires that we subtract from the total energy, the total 
mass of the system at infinite separation (i.e. the sum of 
the masses of the individual black holes) . Unfortunately, 
there is no rigorous way of defining the total mass of an 
individual black hole at finite separation. This is a fun- 
damental uncertainty in defining a relativistic effective 
potential. 

The total mass of an individual black hole includes 
both its irreducible mass and the kinetic energy associ- 
ated with the spin. In general, we cannot rigorously de- 
fine the spin of an individual black hole, nor do we know 
how to compute the total mass of an individual black hole 
given the spin and apparent-horizon mass. A reasonable 
approximation is to define the total mass of an individual 



black hole, M, via the Christodoulou formula[22 

52 
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(7) 



where S is the magnitude of the spin of the black hole 
and we can approximate the irreducible mass, M^,,, with 
the apparent-horizon mass. 

To avoid some of these ambiguities, we will define the 
binding energy as the ADM energy minus the sum of the 
apparent horizon masses. 
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(8) 



Because we keep the apparent-horizon masses constant 
along an EP curve, this constant shift does not affect the 
effective potential. But, since we do not keep the total 
mass constant along an EP curve, the uncertainty in the 
definition of the total mass contributes to the uncertainty 
in the definition of an EP curve. 

Another fundamental uncertainty in our numerical EP 
curves is the assumption that they are constructed of 
turning points in eccentric orbits. When the EP method 
was introduced, it was applied while using the conformal 
imaging decomposition which allowed one to directly set 
the radial momentum of the two holes to zero [6]. In the 
CTS decomposition, there is no mechanism to directly set 
such data. However, as discussed at the end of Sec. [TTl 
we conjecture that our initial data does have zero radial 
momentum. The reasoning behind this conjecture is that 
the time derivative of the conformal metric is still taken 
to vanish, so that at least instantaneously we expect each 
black hole to be at a turning point. 

In order to test the validity of this conjecture, we com- 
pare our numerical effective potentials to post Newtonian 
data. Mora and Will introduced third order conser- 
vative post Newtonian equations for the energy and an- 
gular momentum of a system in terms of an eccentricity e 
and inverse semilatus rectum C. The authors treated the 
black holes as having zero spin and ignored dissipative 
terms. One can find the orbital angular velocity at ei- 
ther pericenter or apocenter through the author's choice 
of definition of e and ^ used in creating these equations. 
In order to include the spin of the black holes, it be- 
comes necessary to include correction terms that create 
equations of the following form. 

Eie,C,u;) = ^AOM(e,C)+£^.eif(e,C,c^) (9) 
+ i;„.co„(e,C,^) + Sspi„(e,Cw) 

J(e,C,^) = AnM{e,0 + S{eX,oj) (10) 
+ JN,co„(e, C, t^) + -'"spi„(e, C, t^) 

In these equations uj represents the spin angular veloci- 
ties of the black holes^ . The self energy and spin terms 



^ It would be preferable to hold the total mass of the individual 
black holes constant. However, as we will discuss in the main 
text, this is not possible. 



^ We use a single parameter lj for simplicity. In general, specifying 
the spin of two black holes would require six parameters. 
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FIG. 2: EP curves for non-spinning equal-mass black holes 
from both numerical and 3PN data. Numerical EP curves are 
plotted as short-dashed(blue, green, and red) lines. 3PN EP 
curves are plotted as solid(black) lines. The Komar sequence 
through the numerical data is plotted as a dash-dot (purple) 
line. The EP sequence through the numerical data is plotted 
as a long-dash(black) line. The EP sequence through the 
3PN data is plotted as a light-solid(orange) line. A boundary 
of the allowable region for the 3PN equations is shown as a 
light-solid(brown) line. 



(-Bocif and S) are derived as expansions of the Kerr for- 
mulae relating mass, spin, and rotational angular veloc- 
ity. The Newtonian correction terms (-En.coh and Jn.coh ) 
stem from the conversion of total mass to irreducible 
mass and the "Spin" terms (i^spin and Jspi„) represent 
spin-orbit effects. The parameter space of these equa- 
tions are shown below in Figs. [S] and [51 

With this parametrization, one can easily construct 
3PN EP curves using straight-forward one-dimensional 
root finding. To create a sequence of constant angu- 
lar momentum, one can apply a sequence of values for 
the eccentricity to Eq. (jlOp. and for each find the cor- 
responding C that returns the desired value of angular 
momentum. One can then use those values of e and C 
in Eq. ^ to find the energy of the system. In Fig. [21 
we plot for comparison the EP curves for non-spinning, 
equal-mass binaries from both the numerical data and 
the 3PN equations. The energy is plotted as a function 
of the dimensionless orbital angular velocity mfl, where 
small mfl corresponds to large orbital separation. In- 
cluded on the graph are data from the Komar sequence 
and the EP sequences extracted from the minima of the 
3PN EP curves and the numerical EP curves. 

At large values of angular momentum (large orbital 
separation), the solid 3PN data agrees well with the 



dashed numerical data. The numerical data and the 
3PN data begin to diverge as the angular momentum 
decreases. This is not surprising as it is well known the 
PN expansion is less accurate for tighter binary systems 
and the same is true for the numerical initial data mod- 
els. It seems from the good agreement between numerical 
and 3PN EP curves that the numerical data we construct 
using the CTS approach are reasonably close to turning 
points, and seem to asymptote to turning points as the 
system becomes more Newtonian. 



IV. MEASURING ECCENTRICITY 

The comparison of our numerical EP curves with 3PN 
EP curves in Sec. IIIII gives us confidence that the numer- 
ical data represent systems at either pericenter or apo- 
center. Since each numerical EP curve is constructed 
from a sequence of models where the mass, angular mo- 
mentum and spins are held constant and are at turning 
points, it is reasonable to use these curves to try to de- 
fine the eccentricity for a given model. For Newtonian 
binaries, eccentricity can be measured using relative sep- 
aration at pericenter dp and apocenter da'- 



da 



da + dp 



(11) 



We use ed to denote the eccentricity of numerical models 
based purely on separations. 

By applying Newtonian equations of motion, one can 
replace the above definition of eccentricity dependent on 
relative separation with a version that is dependent on 
the orbital angular velocities at pericenter f2p and apo- 
center fla- Not only is the eccentricity found in those 
terms, but the dimensionless inverse semilatus rectum 
can also be found using the same parameters [l^ . 



Co 



c 



^p 4" 



4/3 



(12) 

(13) 



We use en to denote the eccentricity of our numerical 
models as measured using orbital angular velocities. The 
following relationships follow directly from these equa- 
tions: 



c = 



(1 



2/3 



m^a 

(W 



2/3 



(14) 



We note that these equations (at either pericenter or apo- 
center) are the only places in any of our work where the 
sign of the eccentricity matters and so for the remain- 
der of the paper we will ignore the middle (pericenter) 
relationship and assume eccentricity is negative at peri- 
center. This has the added benefit of cleaning up many 
of the figures below. 
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A. Direct Measurement 

Our next goal is to determine the eccentricities of the 
data in our EP curves and determine if the definition 
of eccentricity in Eq. ((T^) yields reasonable results. We 
start by computing cq for the non-spinning equal-mass 
models represented in the numerical EP curves displayed 
in Fig. [2] (also seen in Fig. [1]). The measured eccentric- 
ities are shown in Fig. [31 plotted against mSl. The eo 
definition of eccentricity should give reasonable results 
when the corresponding orbit is sufhciently Newtonian. 
The orbits become more relativistic as the total angu- 
lar momentum associated with the orbit gets smaller. 
Throughout this paper, we will use the following con- 
vention to easily differentiate which EP curves represent 
large, intermediate and small values of angular momen- 
tum. EP curves with dot-dashed lines represent large 
values of angular momentum, solid lines denote the mid- 
dle range of angular momentum, and the dashed lines 
have small values of angular momentum. 



Small J/|a,m 
Intermediate J/|a,m 
Large J/|a,m 




0.3 



FIG. 3: Eccentricity of non-spinning equal-mass black holes 
computed along numerical EP curves. The en definition of ec- 
centricity is plotted against the orbital angular velocity. Neg- 
ative values of en correspond to models at pericenter, while 
positive values correspond to apocenter. Large values of the 
orbital angular momentum are plotted as dot-dashed(blue) 
lines, intermediate values as solid(green) lines, and small val- 
ues as dashed(red) lines. 

Figure [3] shows only a limited range of eccentricities 
for each value of angular momentum. There are sev- 
eral reasons for this, both physical and computational. 
First, to compute the eccentricity, we require data from 
two corresponding turning points on the same EP curve. 
That is, we need two points with the same value of the 
binding energy. Because we do not construct models at 



arbitrarily large separation, some data at pericenter have 
no matching data at apocenter. In this case, the eccen- 
tricity cannot be computed. Clearly, we cannot compute 
eccentricities using this method for pericenter data cor- 
responding to unbound orbits, but there there are addi- 
tional limitations associated with the shape of the effec- 
tive potential at small separation. Because of strong-field 
effects, the effective potentials reach a local maximum at 
small separation (cf. the effective potential for massive 
test particles orbiting Schwarzschild) . For bound orbits, 
we are limited to computing eccentricities for configura- 
tions with energies lower than that of the local maximum 
for each EP curve. So, the eccentricities plotted in Fig. [3] 
correspond to data in the neighborhood of the local min- 
imum of each EP curve and extending only as far as the 
lowest local maximum on either side. 
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FIG. 4: The relative difference between the en and ed defi- 
nitions of eccentricity for the non-spinning equal-mass black 
hole numerical EP curves. Lines as in Fig. |3l 

To test whether the numerical data is behaving as ex- 
pected, we turn to the definition of eccentricity found in 
ed- This definition should be reasonable for large separa- 
tions, but will break down as the coordinate separation 
of the two black holes decreases because of coordinate ef- 
fects near the black holes. In Fig. 21 we show the relative 
difference of the two measurements (en — ed)/en- Note 
that we have manually removed data corresponding to 
points near the minima of the EP curves since both def- 
initions of eccentricity yield zero at the minimum and 
the relative error for neighboring points is dominated 
by numerical noise. However, it is easy to find where 
those points would have been as the different lines be- 
come somewhat jagged in the region of zero eccentricity. 

Figure 31 shows the behavior we would have expected. 
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For large separations (large J / fim) the coordinate separa- 
tion d and the orbital angular velocity should both yield 
reasonable estimates of the eccentricity and we see that 
the relative error is tending to zero as J//im increases. 
Clearly, the gauge dependence of the coordinate separa- 
tion d will cause to become less reliable for smaller 
separations (small J/ ^m) and indeed, we see the relative 
difference increase as J/fj,m decreases. Because the or- 
bital angular velocity fl is gauge independent, we expect 
that eo will yield a better definition of the eccentricity, 
however, we need an independent standard against which 
we can measure the reliability of en- 



B. Post-Newtonian Measurement 
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FIG. 5: Parameter space of J(e, (^) (Eq. (I10|l 'l for non-spinning 
equal-mass black holes. Short-dashed(green or black) lines of 
constant e are constructed using the 3PN equations. The 
solid(red) line represents all 3PN models with a specific value 
of m£l = 0.0614 corresponding to a particular numerical 
model. The horizontal dotted(black) line represents the cor- 
responding orbital angular momentum of that model. 



To test en in the more relativistic regime we return to 
post-Newtonian theory. The 3PN equations for the en- 
ergy and angular momentum in Eqs. ([5]) and (jlOp can be 
used to compute the eccentricity of initial data in sev- 
eral ways, all using information from a single initial-data 
configuration[12, 23, 24]. The two most useful ways are 
based on using values for either E\, and SI, or J and 
from a given initial-data set. Both methods yield similar 
but distinct values for the eccentricity. We will denote ec- 
centricities obtained using the energy via tE^^ and using 
the angular momentum via ej. 
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FIG. 6: Parameter space of E{e,() (Eq. Q) for non-spinning 
equal-mass black holes. Lines as in Fig. [5] except the horizon- 
tal dotted(black) line represents the corresponding binding 
energy of the numerical model. 



To use Eqs. © and (fTU|) to find eccentricity, we need 
to simplify the dependencies. Those equations depend 
on the eccentricity e, the inverse semilatus rectum C, and 
the individual black- hole spins (represented by to). The 
spin dependence can always be fixed. For now, we con- 
sider non-spinning black holes. Next we apply Eq. ([14]) 
to replace e with mfl and C. The equations now depend 
on only mil and To find a 3PN value of C, we set one 
of the equations (say the energy equation) to the con- 
stant (energy) taken from an initial-data set and replace 
mQ with its value from the same data set. We can then 
use one-dimensional root finding to obtain a value for 
Finally, using Eq. ([14]) again, we can obtain e. 

There is an issue when using root-finding methods on 
the modified equations. Eqs. @ and (fTU]l are polyno- 
mials with multiple roots, so we must determine which 
value of ( to use. Figure |S] shows the parameter space 
for the 3PN angular momentum from Eq. (jlOp for the 
case of non-spinning equal- mass black holes. The angu- 
lar momenta for lines of constant eccentricity are plotted 
against C as dashed lines. A solid line shows all 3PN con- 
figurations having a constant value of mfi = 0.061355. 
This value of mSl was chosen because it corresponds to 
the minimum of one of the numerical EP curves. The 
horizontal dotted line displays the angular momentum 
from the numerical data set. As can be seen, there are 
two values of C where these two lines intersect, and both 
correspond to valid roots of the equation. The small- 
est positive root corresponds to a very small eccentricity 
which we would expect for the given data set. The second 
smallest positive root yields an eccentricity somewhere 
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between ej = 0.25 and e,/ = 0.50. It is unlikely that the 
minimum of an EP curve would have such high eccen- 
tricities and negative values of C aren't allowed, so we 
always choose the smallest positive root for The 3PN 
energy from Eq. ([9]) yields similar results as can be seen 
in Fig. \6\ Again, we always choose the smallest positive 
root for (. 



Small J/\im 
Intermediate 




0.2 



FIG. 7: Eccentricity of non-spinning equal-mass black holes 
computed along numerical EP curves. The 3PN e_B^ definition 
of eccentricity is plotted against the orbital angular velocity. 
Lines are as in Fig. |3] The x symbols mark the minima of 
each EP curve. 

In Fig. [7] we show the eccentricity ce^ of the same EP 
curves considered in Fig. [3] but computed using the 3PN 
energy equation. The minima of the EP curves, which 
would have en = 0, are marked with x 's. Qualitatively, 
these results resemble the previous direct measurements 
shown if Fig. [31 however there are differences. First, note 
that the the minima of the EP curves do not correspond 
exactly to e^t =0 (as first noticed in Refs.p^. [23l.l23|). 
Also, for high values of J/fim, can be evaluated for 
more of the data points on these EP curves than is possi- 
ble for eo . The evaluation of was limited in this range 
because the numerical data did not extend out to suffi- 
ciently large separations and eo requires pairs of turning 
points to measure the eccentricity. Because requires 
information from only a single initial-data set, it can be 
computed for some data points where cannot. While 
not presented here, using the angular momentum equa- 
tion to compute ej delivers results that are qualitatively 
similar. 

Finally, our goal has been to gauge whether or not eo 
was a reasonable definition of eccentricity. In Fig. [51 we 
show the difference between the en and 3PN definition 
of eccentricity eE^- We plot en — e^^ against en rather 
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FIG. 8: Difference between en and es^ plotted against en for 
non-spinning equal-mass black holes computed along numer- 
ical EP curves. Lines are as in Fig. [S] 



than a relative difference to avoid division by small num- 
ber issues that make the graph difficult to read (recall 
that en and do not evaluate to zero for the same 
data points). As expected, there is better agreement for 
more Newtonian configurations (large J/ fim) which di- 
minishes as J//im decreases. There is some jaggedness 
at en = 0.0 caused by the polynomial fitting used to 
estimate the minima of the EP curves. A careful exam- 
ination of the apocenter data points (positive en) shows 
that relative errors for modest values of en do not exceed 
20% for even the smallest values of J//im. Of course, the 
relative errors near en = are unbounded. This com- 
parison suggests that the use of either en or eEf, yield 
reasonable measures for the eccentricity when applied to 
the numerical initial-data sets, although we should be 
more cautious in trusting results for pericentcr data with 
small values of the angular momentum. 



C. The Komar-Mass Difference 

The first application of the ee,, and ej definitions to 
non-spinning, equal-mass black-hole binary initial data 
was undertaken by Berti et. al.jl^. In this work, they 
considered initial data sets that satisfy the Komar-mass 
criteria for circular orbits and showed that the 3PN def- 
initions of eccentricity ce^ and ej yield non-zero results 
for these configurations which are supposed to be in cir- 
cular orbits. We reproduce these results in Fig. [51 We 
note that, while it is true that the circular orbit data 
have non-vanishing 3PN eccentricity, the magnitude of 
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this eccentricity smoothly approaches zero as the separa- 
tion increases. We also note, as pointed out in Refs. [12 
and [2^, that the energy based definition e^b yields con- 
sistently smaller values of eccentricity for the "quasi- 
circular" data than does the angular momentum based 
definition ej. 



O.li 




FIG. 9: The 3PN eccentricity measures es^ and ej applied 
to both Komar and EP sequences of non-spinning equal-mass 
black holes. Note that the EP minima yield quasi-circular 
data with a smaller eccentricity than is obtained from the 
Komar-mass ansatz. 

In Ref. we showed that circular orbits defined 

by the EP method yield models that are very similar 
to those defined by the Komar-mass ansatz. We also 
showed that, while very similar, all quantitative mea- 
surements of the quality of the circular orbits showed 
that the EP method yields better results. This is again 
true if we compare the 3PN eccentricities computed for 
non-spinning, equal-mass binaries in circular orbits de- 
fined by the EP method and the Komar-mass ansatz. 
The results are shown in Fig. [51 where it is clear that 
the the EP method yields consistently smaller values of 
eccentricity for quasi-circular orbit data. 

Although the EP method yields consistently better 
results for circular orbits than can be obtained using 
the Komar-mass method, the differences are in general 
not sufficiently significant to outweigh the considerable 
computational expense associated with the EP method. 
When we compute eccentricities using the eo definition, 
the overhead of using the EP method is even larger. It 
would be useful to find another means of estimating the 
eccentricity of binary initial data. 

In Fig. [ini we again plot the eccentricity en of the same 
set of non-spinning, equal-mass initial data. However, 
on the horizontal axis, we plot the dimensionless Komar- 




FIG. 10: The eccentricity measure en applied to non-spinning 
equal-mass black holes computed along numerical EP curves 
and plotted against the Komar-mass difference AK. Lines 
are as in Fig. [3] 

mass difference 

AK = {Eadm ~ Mk) (15) 

We see a very strong correlation in the data, though 
the correlation weakens as we move further from quasi- 
circular orbits. It is worth noting that the EP curves 
with low J I puTfi curve back towards zero AK for con- 
figurations with negative eccentricities (pericenter) and 
small J I iim. This behavior is for data in the region of 
the local maxima in the EP curves. That this occurs is 
consistent with the notion that these local maxima rep- 
resent unstable quasi-circular orbits, and hence we will 
find the Komar-mass difference vanishing in this region. 
More importantly, it shows that all of the definitions of 
eccentricity we have used will break down in this highly 
relativistic region. 

V. COROTATION 

So far we have restricted ourselves to the case of non- 
spinning black holes. However, corotating configurations 
have received considerable attention in spite of the fact 
that we do not expect to see corotating black holes in 
nature. In addition to the non-spinning case, Berti et. 
al. [23I also computed the 3PN eccentricities for the coro- 
tating data presented in Ref. [Toj and we reproduce these 
results in Fig. [TT] 

One of the primary reasons that the initial data for 
corotating black holes has been studied so extensively 
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is that there is a simple and unique means of enforcing 
the condition of corotation on the black holes. This is 
in contrast to any attempt to enforce a specific value of 
spin (even no spin) on each black hole, which necessarily 
includes the uncertainty in how we define the spin of an 
individual black hole in a binary configuration. 
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FIG. 11: The 3PN eccentricity measures es^ and ej applied 
to both Komar and EP sequences of corotating equal-mass 
black holes. The dot-dashed(blue) lines correspond to quasi- 
circular data defined via the Komar-mass ansatz and where 
the Newtonian notion of corotation is used in the 3PN equa- 
tions. The dashed (red) lines show the improvement obtained 
by including the IPN correction to the notion of corotation. 
The solid(black) lines show the added improvement of using 
quasi-circular data based on the EP method. 

The Newtonian concept of corotation implies that each 
black hole rotates with a spin angular velocity ujs that is 
equal to the orbital angular velocity. This Newtonian 
notion of corotation (i.e. ujs = f2o) was used to fix the 
spin parameters in Eqs. ^ and (fTO|) when using these 
3PN equations to compute es^ and ej. To our knowledge, 
this Newtonian notion of corotation has been used in all 
PN computations dealing with corotation (cf Refs. [H, 
[IS, [23, 111). However, in Ref. 0, we have shown 
that there are relativistic corrections to the spin angular 
velocity associated with corotating black holes. We find 
Ws, including the IPN correction, to be of the form 



= r!o(l-^(mr!o)'/' + •••), 



(16) 



where 77 = /i/m is the symmetric mass ratio which takes 
the value of 1/4 for equal-mass binaries. Also shown in 
Fig. [11] are the two 3PN eccentricities computed using 
the corrected definition for cus- We find that using this 
corrected definition significantly decreases the 3PN esti- 
mated eccentricity for these circular-orbit models. 



As with the non-spinning case, we can also compute 
the eccentricities for the corotating equal-mass binaries 
in quasi-circular orbits defined in terms of the EP method 
rather than the Komar-mass ansatz. Including also the 
improved definition for uJs in the 3PN definitions of ec- 
centricity, we find the evaluated eccentricities are small- 
est when evaluated for quasi-circular data defined by the 
EP method. This can also be seen in Fig. [TTl 

Even with the correction to cj^, we notice the magni- 
tude of the 3PN eccentricities computed for "circular" 
data are consistently larger for corotating binary data 
than for non-spinning data. We cannot be certain why 
this is the case, but we should keep in mind that there 
is an inherent inconsistency in any attempt to attach the 
notion of eccentricity to corotating configurations. All of 
our definitions for eccentricity, including the 3PN defi- 
nitions that can be evaluated using information from a 
single data set, ultimately rely on information from both 
a pericenter and an apocenter configuration. For eccen- 
tric orbits, the spins of corotating black holes will change 
from pericenter to apocenter. Since the spins (includ- 
ing spin-orbit and spin-spin couplings) contribute to the 
total energy, this variation of the spin throughout the 
orbit must impact upon our definitions of eccentricity. 
This is likely to cause few problems when considering 
nearly circular orbits (where the spin varies little from 
pericenter to apocenter), but our definitions of eccen- 
tricity may not yield reasonable results for orbits that 
deviate significantly from being circular. 



VI. CONCLUSIONS 

In this paper, we have examined several basic ques- 
tions associated with the construction of binary black- 
hole initial data. When we construct binary initial-data 
sets using the extended CTS equations and fix the freely- 
specifiable parts of the data and boundary conditions 
to be consistent with the assumptions of quasiequilib- 
rium, we obtain models for black holes in circular orbits. 
But what happens if we set aside the quasiequilibrium 
assumption that imposes circular orbits? The resulting 
data can no longer evolve in a quasiequilibrium manor as 
the orbit will have a significant eccentricity. Our inves- 
tigations suggest that the initial data models we obtain 
represent, in general, binaries that are at turning points 
(either apocenter or pericenter) of some general eccentric 
orbit. 

The specific notion of quasiequilibrium that we set 
aside is implemented by imposing either the Komar-mass 
condition or by choosing the minimum of an EP curve as 
our circular orbit model. If, as it seems, general initial- 
data models on an EP curve are at turning points, then 
we can use information from these models to estimate 
the eccentricity of the model's orbits. Of course, there is 
no unique definition of eccentricity. We have examined 
several possible definitions for eccentricity. Using only 
information from the initial-data sets on an EP curve, 
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we have defined two eccentricities for an orbit (eo and 
e^), but these definitions require that we have representa- 
tive models at both the pericenter and apocenter turning 
points of a given orbit (assumed to have constant binding 
energy). We have compared these definitions of eccentric- 
ity to the 3PN definitions (eg^ and ej) developed by Mora 
and Willfl^. All of the definitions agree quite well for 
non-relativistic orbits. They are also in reasonably good 
agreement for more tightly-bound and relativistic orbits 
as well, although the results are quantitatively different. 

One might ask which definition is better for more rel- 
ativistic situations. However, it isn't clear that there is a 
meaningful answer. In comparing the numerical and 3PN 
eccentricity measures for relativistic cases, we have no- 
ticed an interesting and unexpected feature of the 3PN 
equations. Figure [5] compares numerical and 3PN EP 
curves for equal-mass non-spinning black-hole binaries. 
The EP curves each have a constant value for the orbital 
angular momentum. All of the numerical EP curves cover 
a finite range of separations (parameterized by mflo). 
This is because it becomes computationally expensive to 
compute models at very large separations (small milo) 
and it becomes increasingly difficult to obtain conver- 
gent solutions at very small separations. What was un- 
expected is that the 3PN EP curves also cover a finite 
range of separations. 

For sufficiently small values of J/nm, the 3PN EP 
curves do not extend to configurations with arbitrarily 
large separation. To make this clear, Fig. [5] includes a 
curve that marks the boundary (for both large and small 
separation) of the 3PN EP curves. This is most easily 
seen for the 3PN EP curves near the bottom of Fig. [2l 
Here we can see that the 3PN EP curves do not extend 
to arbitrarily small values of mQo , but the numerical EP 
curves do (although we do not compute them for arbi- 
trarily large separation). 

That this behavior is not an artifact of our method 
for computing the 3PN EP curves can be seen by ex- 
amining Fig. [51 Recall that angular momentum is held 
constant along EP curves, so an EP curve is repre- 
sented by a horizontal line in this figure. Notice that 
all curves of constant eccentricity have a local minimum. 
For J/ fim > 3.4, an EP curve can extend from e = ^ 1. 
However, for J/ iim ^ 3.4, can only extend from e = to 
a maximum eccentricity that is less than 1, and which is 



determined by which curve of constant e has its minimum 
tangent to the EP curve. 

It is physically reasonable that large eccentricity or- 
bits should not exist as the orbits become sufficiently 
relativistic since the finite size of the black holes would 
lead to a collision. However, the PN equations treat the 
black holes as point particles and should not be sensitive 
to this limitation. Also, we remind the reader that the 
numerical EP curves have no difficulty in extending to 
large separation (e 1) for relativistic orbits. This sug- 
gests that, for relativistic orbits, there may be a problem 
with the definition of eccentricity and semilatus rectum 
used by Mora and Will[l3| to parameterize the PN en- 
ergy and angular momentum. We do not consider this 
to be a serious problem, but rather an indication that 
any definition of eccentricity is of limited value for orbits 
where the interaction between the black holes becomes 
highly relativistic. 

Given the strong correlation between the various mea- 
sures of eccentricity and the difference between the Ko- 
mar mass and the ADM energy as measured by AK and 
seen in Fig. [TOl we suggest that perhaps AK can serve 
as a useful invariant means of parameterizing the eccen- 
tricity of an orbit. However, initial-data studies can at 
best suggest possible useful parameterizations. It will be 
most useful to evolve initial data that are significantly 
eccentric and examine the orbital dynamics to better un- 
derstand both the parameterization of eccentricity and 
its effects on the dynamics. For example, it may be par- 
ticularly interesting to explore the evolution of eccentric 
initial data for an "orbit" that has no pericenter turning 
point on its EP curve. 
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